Long Range Correlations in the Disordered Phase of a Simple Three State Lattice Gas 
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D , We investigate the dynamics of a three-state stochastic lattice gas, consisting of holes and two 

■ oppositely "charged" species of particles, under the influence of an "electric" field, at zero total 
charge. Interacting only through an excluded volume constraint, particles can hop to nearest neigh- 
bour empty sites. With increasing density and drive, the system orders into a charge-segregated 
state. Using a combination of Langevin equations and Monte Carlo simulations, we study the 
steady-state structure factors in the disordered phase where homogeneous configurations are stable 

■ against small harmonic perturbations. They show a discontinuity singularity at the origin which in 
O ' real space leads to an intricate crossover between power laws of different kinds. 

PACS: 64.60Cn, 66.30Hs, 82.20Mj 

Keywords: Driven lattice gas; Monte Carlo simulations; Langevin equations; Long-range correla- 
te , tions. 
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. For a system with short range microscopies, placed in thermal equilibrium, long range spatial correlations are 
absent in general. Their presence is a typical signal that the system is at a singular point. On the other hand, for 
r : systems in non-equilibrium steady states, long range correlations are often observed |l[]. With conserved dynamics 
' ^ ' they generically arise if the system has strong spatial anisotropics and lacks detailed balance. A good example is the 
driven Ising lattice gas |^,|| . The role of the drive ( "external field" ) is to bias hopping rates along a particular direction 
on the lattice. In addition to many other unexpected features such as non-Hamiltonian fixed points controlling the 
critical behavior jdj , the system exhibits long range spatial correlations at all temperatures above criticality, as a result 
of the breakdown of the fluctuation-dissipation theorem (FDT). In momentum space, this appears as a discontinuity 

■ singularity of the structure factor at the origin ||] . 

, One can generalize the model of ref |^ in a fashion similar to the one leading to the Ising model to spin-1 |Q or 
' Potts models, by considering more than one species of particles. To start with, we have two species of particles 
(referred as -l-'s and — 's) driven in opposite directions and holes (empty sites). To keep the model simple, however, 
' we may neglect the usual Ising nearest neighbor interaction and keep "only" the excluded volume constraint and the 
bias. This multi-species model, in both one and two dimensions, has been thoroughly investigated l^-p^j. In the 
simplest scenario charge (H — ) exchange is not allowed. Monte Carlo simulations [^j in two dimensions and mean-field 
studies ||,|l^ show that there is a transition, controlled by particle density and drive, from a spatially homogeneous 
(disordered) phase to a charge segregated one, where the excluded volume constraint leads to the mutual blocking of 
particles. If we soften the excluded volume constraint - by allowing exchange between nearest neighbor, oppositely 
charged particles on a much slower time scale than the particle hole exchange - the above transition still survives . 

■ In the disordered phase, these two versions behave in much the same way. Thus, for simplicity, we will study only the 
Ch ! case without charge exchange in this paper and focus on the spatial correlations and structure factors. 
O ' We consider a two dimensional fully periodic lattice with L± x L|| sites, each of which can be empty or occupied by a 

single particle. Since we have two species we need two occupation variables n+ and n~, with n being or 1, depending 
on whether a positive or negative particle is present at site x. Although we refer to these particles as "charged" , they 
do not interact via the Coulomb potential. Instead, they respond to an external, uniform electric field E, directed 
along the -t-siy-axis. We restrict ourselves to zero total charge, i.e. X^xl^x ~ '^x] = 0. In the absence of the drive 
' the dynamics does not distinguish between the different species: both types hop randomly to nearest neighbor empty 
sites, with a rate P. The electric field introduces a bias into these rates in such a way that jumps against the force 
will be exponentially suppressed. During one Monte Carlo step 2L_lL|| nearest neighbor bonds are picked randomly. 
If a particle-hole pair is encountered, an exchange takes place with probability W = P min{l, exp(g£' (5a;|| )}, where 
g = ±1 is the charge of the particle and (Sxy = ±1 is the change of the coordinate of the particle due to the jump. 

For our simulations, we set P = 1. Using 30 x 30 and 60 x 60 lattices, we initialized the system with random 
configurations of various particle densities and carried out runs ranging from 2.5 to 5 x 10^ MCS. After allowing 
62500 MCS for the system to settle into a steady state, we measured the Fourier transform of every 125 MCS. 
The structure factors are then constructed from averages of these measurements, according to 
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5"^ (k) = -^^e-^ (n>^), (1) 

where a,/? = +, — and k = (^X7^'~TJ^) ^ ■ Note that S"^^ is the Fourier transform of the usual equal-time 

correlation function G"''(x) = {n'^n^) — {n'^){n^)- Thus, if G is even in x, S will be real. Equivalently, an imaginary 
part of 5 corresponds to the part of G which is odd in x. Due to charge symmetry, we expect G'^^ = G . Clearly, 
both must be even in x, so that the associated 5's are real. On other hand, due to the drive, we have 

G+-{^,E)^G+~{x^,-x\\ ,-E) (2) 

so that S'^ may have an imaginary part (which must be odd in E). Finally, G ^(x) — G^ (^x) is just a mathe- 
matical identity. 

In Fig. ^ we present the results for the three independent 5's found in the larger system. Before discussing the 
data in detail, we will first present the theoretical framework within which they can be understood. In particular, we 
will illustrate the emergence of discontinuity singularities in the structure factors at k = 0, and their consequences for 
long-range correlations in real space. This will then be followed by a comparison between our theoretical predictions 
and the simulation results. 

Since our interest lies in the behavior at large distances (or small k), the simplest approach relies on continuum 
field theory. Starting from the dynamics at the microscopic level, specified by the above rates and the associated 
Master equation, there are several ways to arrive at a coarse-grained description, in terms of equations of motions 
for the local i-dependent densities. The most systematic one is to perform an Sl-expansion pj], where one splits the 
particle densities describing a block b (size O), around point x, into a macroscopic part (gi^and a fluctuating one 



a ^ "± =£<±(x,i)-ff]-V2^±(x,i). (3) 
x' eb(x) 

For generality, we consider the d-dimensional case when X|| is directed along the electric field and x^ is in the d-1 
dimensional subspace, perpendicular to the field. After taking the continuum limit the result for the macroscopic part 
is a set of mean-field equations of motions for the slowly varying local densities: 

dtQ^ = -vr{e± V (1 - - ^?-) ±ex||^?±(l - - q-)Y (4) 

where 



r 

r,i 



(5) 



II 

is the diffusion-matrix. is diagonal and isotropic in the d-1 dimensional subspace, thus characterized by a number 

r_L. V is the asymmetric gradient operator, namely for any two functions / and g, / V g = /Vg — gV/. e = 
2tanh(£'/2) is the coarse-grained bias and X|| is the unit vector along the a;|| direction. 

Since (^) is a continuity equation, it trivially admits t-independent solutions which are homogeneous in space: 
0^ (x, t) = g . For simplicity, and to ease comparison with simulation data, we have chosen equal densities for 
both species. At the mean-field level, this solution describes the steady state. The correlations of interest are then 
associated with x^, the scaled deviations (^) from the mean-field steady state. 

Returning to the il-expansion, the result for the fluctuating part, at lowest order, is a Fokker-Planck equation. For 
our purposes, the equivalent Langevin equation is more transparent. At this order, its deterministic part is linear and 
the (conserved) noise is Gaussian. Focusing on the fluctuations about the homogeneous phase, the result is: 

^*U-(x,i))-^(^Mx-(x,i))-i V^-(x,t)) ' (6) 



where 



_(C++{V) c+-{v)\_ 

(1 - ^)vrv - (1 - 3^)er||a|| ^vrv + geT^^dw 

gVrV - geT^id^i (1 - ^)VrV + (1 - 3^)er||5|| 



(7) 
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Here, ri^{x,t) are Gaussian, white noise terms, satisfying: 

(77r(x,i)) = , «(x,t),7f(x',i')) = 2S^^a,,S{^-^')S{t~t') (8) 
a, = ; i,j = 1,2,.. .d . 

{<7ij) = cT is the noise matrix: 

Similar to Fj^, crj^ is diagonal and isotropic in the d-1 dimensional subspace, characterized by a number a^. In 
equilibrium systems, the FDT guarantees cr cx; F. In particular, in the absence of the drive, we would have cr = 
p{\ — 2p)F here. However, when driven, the proportionality is not expected to hold, in that the diffusion and noise 
matrices would be renormalized differently by the drive e. This is certainly the situation in the driven single species 
case [Q. Finally, we point out that 77+ and ri~ are uncorrelated, due to the fact that charge exchange is not allowed. 

To find the correlations and structure factors from eqns. (^^), we introduce the Fourier components for the 
fluctuations: 

X^{K^) = j dtd'^x X^(x,t) e-*("*+'^'') , (10) 



and similar ones for the noise, so that 

a,(3 = +,- ; i,j = l,2,...d 
Then the solution to (o) is trivial: 



ivTiK^)) - , (ryf(k,..)r;f (k',..')> = 2<S"'3fT,, [(27r)'^+i<5(k + k')<5(c. + c.')] (H) 



^tJt'^n = iC{^)~^u;)-' ( 'i'^-^M] . (12) 

Note that, in k space, (£++,£^^) and (>C^ ,C ^) are complex conjugate pairs. 

As expected, (x^(k, cj)) = 0, since these are the fluctuations about the conserved average densities. Their correla- 
tions are just the dynamic structure factors, easily obtained from (fO) and (p^. From the definition: 



(X"(k, c.)x''(k', Lo')) ^ S'^'^ik, Lo) [{27rr+'d{k + k')Siuj + co')] (13) 
we find the two independent iS"s: 

2kcrk 



5++(k,c.) = yj^^^^ {\^^- \' + I £+- 1^} 

S+-{k,u;) = ^^^{-2C--{^k)C+-{^k)} (14) 
I A(k,u;) 1^ 

where 

A(k, uj) = det{C{ik) ~ iuj) . (15) 

To find the equal-time correlations, we need the steady state structure factors, which can be obtained from (jlj) by 
an integration over lo . The results are: 

_ kak I C++{l^^) |2 



-iTr/:(ik) detC{ik) 
S+-(k) -C-i^)C+-i^k) 

* ^^"-iTr/:(zk) det£(zk) ' ^^^^ 

Note that —^TTC{ik) = (1 — p)kFk is positive definite, while we need to require det £(?k) > for all k 7^ to assure 
that the system is within the linear stability boundary. Finally , using (|^), eqs. (|l6|) take the explicit form: 
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5++(k) 

Re{5'+-(k)} 
Im{5+-(k)} 



(1-^) kak (krk) 



(l-3g)'^ p2^2L2 



(1 - 2^) krk (krk)2 + (1 - 4e)r|e2fc2 
^ fecrfc (krk)2 — ^i^^-j r|e^fc| 



(1 - 2^) krk (krk)2 + (1 - Ag)rje'^k1 
2g fc(Tfcr||efc|| 



(17) 



(1 - g) (krk)2 + (1 - 4e)r2e2A:2 



A key feature of these structure factors is that, unhke in equihbrium cases, ah are singular at the origin. For 
both 5++ and Re{5"' }, this singularity is exhibited as a discontinuity, i.e., limfcu^o S{0, k\\) ^ limkj^^o 5'(kj_, 0). In 
particular. 



and 



limfc„^o 5++(0,fc||) 
linik^^o5++(ki,0) 

limfc„^oRc{g+-(0,fc||)} 
limk^^oRe{5+-(ki,0)} 



ril cr_L 



(l-3e-)2 



(1-^)2(1-4^) 



rii CT± 



(1 - m 



(l-^)(l-4^) 



(18) 



(19) 



In general these ratios are not unity. Note that these singularities come not only from the generic FDT-brcaking 
property, namely ^ 7^ but also from the specifics of this particular driven system. This second factor is a 

monotonically increasing function of g reaching 00 at ^ =1/4. This divergence is related to the instability of the 
homogeneous phase, as we will sec below. 

On the other hand, though ImjS'^ (k)} vanishes for k — > in any direction, discontinuities are present in higher 
derivatives. In configuration space, these singularities translate into power law decays of the equal-time correlation 
functions (x"(x' + t)x''(x', f)), which is, thanks to translational invariance, independent of x'. More precisely, it is 
just 



(2^) 



To carry out the transform, it is convenient define a "mass" m via 

4m2 = (1 -4e)r||£2 



(20) 



(21) 



To keep the system in the homogeneous phase, it is sufficient to impose m'^ > 0. We should note that, in the limit 
eLu — »• 00, the mean-field phase boundary is given precisely by m2 = . Otherwise, for finite sLu, the system does 



not reach the stability limit until g equals i (1 + [27r/£L||]^). Further, it is convenient to rescale the lengths 



X_L 



X_L 

r? 



(22) 



so that F becomes the unit matrix. In terms of these rescaled x, let us define r 
elements of the noise matrix, to keep the notation simple: 



cr_L 



cr± 



|x|. Also, we should rescale the 

(23) 



Now, the transform can be carried out exactly. The results are: 

1 



G++(x) = Vo-V 
G^-(x)=V<7V 



-A 



+- 



1 



rd-2 



( cosh(ma:||) ^— j ^ Kd-2 (r 



G+~(x) = V<tV < B+~sgn(£) sinh(ma;||) ^— j ~^ Kd-2 (mr) 



(24) 
(25) 
(26) 
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where are the parts of G"*" even or odd in being the transforms of the real and imaginary parts of . 

The full correlation is, of course, 

G+-(x)=G^-(x)+G+-(x). (27) 



This decomposition simply reflects the symmetries of the system in the presence of the field (g). In eqns. (|2J-g6|), 
K,j{z) is the modified Bessel function and, provided e ^ 0, the A's and B's are positive constants depending on Fj^ || 
and g. 

The first terms in ( p^ ) and (p5|), bein g p roportional to r~'*, are the well known power law decays ]l5[ | due to 
"FDT-violation" The B-terms in (p4^6|) produce an exponential decay for large r, except along the field. This 
can be easily deduced, since, for large | |, the exponentials of the hyperbolic and Bessel functions cancel. Without 
presenting the detailed asymptotic expansions of (|2^p6|), we illustrate the main results as mr — > oo. 

For x_L / : 



G++(x) , G+-(x) oc ^ 



x^^ - (d- 



(28) 



where the ... represent exponential, short-ranged tails. In this form, we emphasize the three key ingredients of the 
"FDT-violating" power law decays, namely, the dependence on cr|| ^ a± (or cr||/r|| ^ <7±/T±), the dipole amplitude 
and, of course, the r^"^. Also, note that the part of G^ odd in x\\ is seen to be only short-ranged. 

The more interesting limit is for xx — 0. Then, in addition to the above power laws, we have another power: 
I x\\ For all d > I, this power will dominate over the "FDT-violating" component. In fact, if we study d > 3 

systems, even the next leading term in the asymptotic expansion of will be more important than r^''. For typical 
simulations in d = 2, we would write: 

G++(0,X||) , G+-(0,:c||) oc - I +0 (| x,, T^) , 

G+-(0,a;||)cx-sgn(£a:||) l^ii ri -t-O (| x,, T^) . (29) 

Note that, to emphasize the sign of the amplitude of the leading power, we included explicit factors of (—1), so that the 
proportionality constants in (|29| ) are positive. In conclusion, the spatial correlations arc dominated by the expected 
r~'^ power law, except along the field, where a novel | x\\ decay takes over. This new power law comes from 

the coupling between the two species as a result of the excluded volume constraint and the opposite bias. 

Finally, let us turn to comparisons with simulation results. Typically, we find that power law tails are difficult to 
observe, in relatively small systems such as ours. Thus, we focus on the structure factors. Deferring details to a later 
publication, we point out some general characteristics. On the whole, the agreement between data (Fig. ^ and ( p^ ) 
is quite impressive, considering that the latter is only the "first approximation" . For example, 

• in all three cases, S{kj_,0) is independent of k±] 

• the value where Re{5^ (0, fc||)} vanishes agrees well with the theoretically predicted fc|| = y/ {1 — 3p) /{I — p) e; 

• similarly, the ratio Re{S'+^(fc^, 0)}/5++(fc^, 0) is found to be quite close to the expected value of — p/(l — p). 

On the other hand, the data and do not agree so well for, e.g., Re{S'+"(0, y^)}/S++{0, fj). We suspect 
the origin of this discrepancy to be the following. Given that our data were collected quite close to a second order 
transition (so as to observe relatively large fiuctuations), we should expect the longitudinal component of the two- 
point function to suffer considerable renormalization. Certainly, this is the case in the proto model . It would be 
very interesting to take into account finite size corrections, to compute these effects in field theory, and to carry out 
extensive runs for detailed comparison. 

In summary, we have examined the spatial correlations and the structure factors in a simple model of biased diffusion 
of two species. Using both simulation and analytic techniques, we find the expected power law decay, r"'', typical 
of non-equilibrium steady states of a system with anisotropy and subjected to a conservation law. In addition, a 
novel power, r^id-+^)/^^ jg found for correlations along the bias. The general agreement between simulations and a 
simple mean-field theory is surprisingly good, while we await a better theory in order to make detailed quantitative 
comparisons. 
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FIG. 1. Steady state stmeture factors (a) S'++(k) , (b) Re{S'+"(k)} 
and Q = 0.175. Structure factors are plotted against the integer m\\ = -| 
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(c) Im{S+"(k)} for a 60 x 60 system at S = 0.471 



-, while m± 



is taken as a parameter. 
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